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Abstract 

For systems with charged chiral fermions, the imbalance of chirality in the presence of magnetic 
field generates an electric current - this is the Chiral Magnetic Effect (CME). We study the dy¬ 
namical real-time evolution of electromagnetic fields coupled by the anomaly to the chiral charge 
density and the CME current by solving the Maxwell-Chern-Simons equations. We find that the 
CME induces the inverse cascade of magnetic helicity towards the large distances, and that at late 
times this cascade becomes self-similar, with universal exponents. We also find that in terms of 
gauge field topology the inverse cascade represents the transition from linked electric and magnetic 
fields (Hopfions) to the knotted configuration of magnetic field (Chandrasekhar-Kendall states). 
The magnetic reconnections are accompanied by the pulses of the CME current directed along 
the magnetic field lines. We devise an experimental signature of these phenomena in heavy ion 
collisions, and speculate about implications for condensed matter systems. 
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I. INTRODUCTION 


The anomaly-induced transport of charge in systems with chiral fermions has attracted 
a significant interest recently. This interest stems from the possibility to study a new kind 
of a macroscopic quantum dynamics. While the macroscopic manifestations of quantum 
mechanics are well known (for example, superfluids, superconductors and Bose-Einstein 
condensates), so far they have been mostly limited to systems with broken symmetries 
characterized by a local order parameter, e.g. the density of Cooper pairs in superconductors. 
The effects induced by quantum anomalies in systems with chiral fermions are of different 
nature. 

Let us consider as an example the Chiral Magnetic Effect (CME) in systems with charged 
chiral fermions - the generation of electric current in an external magnetic field induced by 
the chirality imbalance [I], see Refs. M for recent reviews and references. In this case, no 
symmetry has to be broken, and the system is in its normal state. However the chirality 
imbalance is linked by the Atiyah-Singer theorem to the non-trivial global topology of the 
gauge field. Since the global topology of the gauge field cannot be determined by a local 
measurement, there is no corresponding local order parameter, and we deal with “topological 
order”. 


This has very interesting implications for the real-time dynamics of a system composed 
by charged chiral fermions and a dynamical electromagnetic field. Indeed, let us initialize 
the system by creating a lump of chirality imbalance localized within a magnetic flux that 


forms a closed loop, see Fig. 1(a) . Magnetic field will induce the CME current flowing along 
the lines of magnetic field B (note that this effect is absent in Maxwell electromagnetism). 
Because the vector CME current acts as a source for the magnetic field, the current flowing 


along B will twist the magnetic flux (see Fig. 1(b) ) and induce a non-zero expectation 
value for the magnetic helicity known since Gauss’s work in XIX century and introduced in 
magnetohydrodynamics by Woltjer [7] and Moffatt [Hj, see also [9J: 


h rn = I cl 3 x A ■ B , 


( 1 . 1 ) 


where A is the vector gauge potential. Magnetic helicity is a topological invariant (Chern- 
Simons three-form) characterizing the global topology of the gauge field. It is mathematically 
related to the knot invariant, and measures the chirality of the knot formed by the lines of 


2 




magnetic field. Because of this, the generation of magnetic helicity will create the chiral 
knot out of the closed loop of magnetic flux - so the topology of magnetic flux will change. 
In this paper we will quantify this statement, and study how the topology of magnetic flux 
changes in real time. We will find that as a consequence of chiral anomaly and the CME, 


the magnetic field evolves to the self-linked Chandrasekhar - Kendall states (see Fig. 1(d) 


). During the evolution, the size of the knot of magnetic flux increases. Moreover, at late 
times this evolution becomes self-similar, and is characterized by universal exponents. 


The evolution of magnetic helicity has been studied previously in the framework of the 
Maxwell-Chern-Simons theory in Refs. PHHT2] (see Ref. ra for study with Maxwell theory). 
The anomaly-driven inverse cascade is discussed in Refs. mm- However the self-similar 
evolution of magnetic helicity has not been reported in these papers. The closest to our 
present study is the paper HU by Tashiro, Vachaspati and Vilenkin , where a simplified 
version of the anomalous magneto-hydrodynamic equation has been applied to cosmic mag¬ 
netic fields. The authors found the power law decay (in terms of conformal time) of the 
chiral chemical potential at the late stage of evolution (see also Ref. [12]). 


We extend the previous studies by elucidating the topology of magnetic flux throughout 
the evolution of magnetic helicity. This is made possible by the use of the eigenfunctions 
of the curl operator in a spherically symmetric domain. Previous studies mm have used 
the eigenfunctions of curl operator in a free-space, i.e. the polarized plane waves. In our 
treatment we can track the magnetic reconnections that transfer helicity from linked to self- 
linked configurations of magnetic flux. We also identify the final state of the system as the 
Chandrasekhar-Kendall state that minimizes the magnetic energy at fixed helicity. 


This paper is organized as follows. In Sec.[TT]we describe the topology of magnetic flux and 
describe the corresponding solutions. In Sec. IH we introduce magnetic helicity spectrum 
and present a qualitative picture of the inverse cascade of magnetic helicity and the role 


played by anomaly. In Sec. IV, we introduce the Maxwell-Chern-Simons equations which 
we will use to study the evolution of magnetic helicity and axial charge density. The results 
of evolution are presented in Sec. [V] In Sec. |VI[ we discuss the relevance of our findings for 
heavy-ion collision experiment. We conclude and discuss possible extensions of the current 
work in Sec. IVIII 
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II. THE CHIRAL ANOMALY AND TOPOLOGY OF MAGNETIC FLUX 


Consider a link /C of N knots of magnetic field with fluxes fa. The corresponding magnetic 


helicity (1.1) of this link is given by 

hm(]G) — fa + 2 <t>i4>j Cij , 


N 


( 2 . 1 ) 


i=l i,j 

where S t is the Calugareanu-White self-linking number, and C tJ is the Gauss linking num- 
bei0 The linking numbers in (2.1) do not always detect the topology of the link; the 
development of the appropriate knot invariants is a very active area of modern mathemat¬ 
ics. The link between the Jones invariant of the knot and Chern-Simons theory has been 
uncovered by Witten m The recent progress includes the HOMFLY knot polynomials, 
Vasiliev invariants, Khovanov and Heegaard-Floer homologies, but the ultimate solution is 


still lacking. In view of this, we will base our discussion on formula (2.1). 


In MagnetoHydroDynamics (MHD), the lines of magnetic field are “frozen” into the fluid, 


and so the magnetic helicity (2.1) is conserved. Moreover, in the absence of dissipation the 


reconnections of magnetic field are absent, and so the topology of the knotted configuration 


is preserved as well - so the two terms in (2.1) are conserved separately. As we will now 


discuss, in fluids with charged chiral fermions the situation changes dramatically due to the 
presence of chiral anomaly. Indeed, the anomaly relation 


d,f A = C A E B (2.2) 

describes the generation of chirality by electric E and magnetic B fields in a topologically 
non-trivial configuration characterized by Chern-Pontryagin number density E ■ B. The 
anomaly coefficient C A for the case of QCD plasma containing N c colors and Nf flavors of 
quarks is given by C A = N c CemG 2 /, with C*em = X)/ qj- The Chern-Pontryagin number 
is easily seen to be equal to the time derivative of magnetic helicity: 

Jt xE . B= -l?£. (2.3) 

This means that when chirality of the fermions is changed, this change is accompanied by 
the change of magnetic helicity, implying the reconnection of magnetic flux. Reconnections 

1 The same formula applies to the helicity of vortex flows, with the substitution of gauge potential A by 
the velocity field v, magnetic field B by vorticity u> = V x v, and the flux fa by the circulation Ki. 
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of magnetic flux in particular can cause transitions between the self-linked (see Figs. |l(c) 


and |l(d) ) and linked configurations of magnetic field described by the first and the second 
terms in ED respectively. Below we will show that such transitions indeed happen as a 
consequence of the anomaly, and the system evolves towards the state in which magnetic 


flux is self-linked, i.e. the entire magnetic helicity is given by the first term of ( 2 . 1 ). 


Before proceeding to the calculations, let us discuss the possible topologies of magnetic 
flux. The Maxwell equations in free space allow for simple solutions with non-zero magnetic 
helicity - these solutions are just circularly polarized plane waves. This is intuitively clear 
since magnetic helicity is parity-odd, and left- and right-circularly polarized waves are the 
simplest P-odd states of electromagnetic field. Since we are interested in describing the 
plasma of a final extent in space, we have however to impose the boundary conditions on 
electromagnetic field. In this case the solutions of Maxwell equations are given by Hopfions 
- configurations in which the loops of magnetic and electric fields are linked. 


On the other hand, the dynamics of electromagnetic fields in the presence of chiral 
anomaly is described by Maxwell-Chern-Simons theory. In particular, the chiral imbal¬ 
ance between the left- and right-handed fermions quantified by the chiral chemical potential 
Ha leads to the generation of electric CME current [I] along the magnetic field: 


Jcme = Ck ha B = a A B , 


(2.4) 


where we introduced the “chiral magnetic conductivity” 04 to allow for the frequency de¬ 
pendence [23]. Unlike the usual Ohmic current, the CME current is topologically protected 
and hence non-dissipative. Hence at late times when Ohmic currents have already dissipated 
away, the r.h.s. of the Maxwell equation V x B = j will contain only the CME current and 
will thus acquire the form 

V x B = a A B . (2.5) 


The solutions of (2.5) have been found independently^] by Chandrasekhar and Kendall 
we will refer to them as CK states, and illustrate their structure in Fig. 1(d) ). It has been 


found by Woltjer [7] that these “force-free” configurations of magnetic field that obey (2.5) 


2 The Editor of their paper [24j wrote: “The results in this paper were derived independently by the two 
authors, and they agreed to write it as one.”. 
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(c) (d) 

FIG. 1. The topology of Abelian magnetic flux: (a) upper left - untwisted loop; (b) upper right - 
twisted magnetic flux; (c) lower left - the self-linked magnetic flux (trefoil knot shown); (d) lower 
right - the self-linked Chandrasekhar-Kendall state. 

minimize the total magnetic energy 


£m = ^ / d x B * 


( 2 . 6 ) 


at a given magnetic helicity (1.1). We thus expect that the CME currents will lead to the 


transition from Hopfion states to CK states at late times, once the Ohmic currents have 
dissipated. We will see below that explicit computations indeed yield this result. 
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III. INVERSE CASCADE OF MAGNETIC HELICITY DRIVEN BY ANOMALY 


A. Magnetic helicity spectrum and Chardrasekhar-Kendal (CK) states 

To discuss the spatial distribution and the inverse cascade of magnetic helicity, let us 
specify the structure of CK states Ref. |23J that are defined as eigen-functions of the curl 


operator, see (2.5). In a free space, a CK state is nothing but a circularly polarized plane 


wave. In this work however, we are interested in EM field in a finite closed system. Let us 
thus consider CK states W^(x; k ) in a spherical domain: 


v x W*(x; k) = ±kW± (*; h) . V • W* (*; h) = 0, 


± 


(3.1) 


where l = 0,1 ,...,m = —l, —l + l,...l. The explicit expressions for W^X] k) in terms 
of spherical harmonics and spherical Bessel functions are given in Appendix. [A] For the 
purposes of the present discussion, we only need to keep in mind that they form a complete 
basis for any divergence-less vector and satisfy the orthogonality relation: 


d 3 x r; k) ■ W^ m ,(x; k') = ^ 5(k - k')5 w 8 mrn ,5 ab , a, b = +, - . (3.2) 


Let us now expand magnetic field B in terms of CK states Wj^ n {x\ k)\ 

^ r 00 r ju 

B(x , t) = J2 / — k2 Wm( k , t)W^(x] k) + ajjk, t)W lm (*; k)] 

l,m ® 

Using V x d = B, we also expand A(x,t) a^J 

^ r°° / u 

A (x,t) = Yl J — k [ a tn( kA ) w L( x 'X) - ai m (k,t) w ^x-,kj\ 


(3.3) 


(3.4) 


Consequently, the magnetic helicity h r 


h rn = I dr x A ■ B , 


(3.5) 


is now given by: 


r°° (ik 

h m (t) = / — kg(k,t ), g(k,t) = g + (k,t) — g~(k,t) 

Jo 71 


(3.6) 


3 We choose the gauge V • A = 0 ; the magnetic helicity h m is gauge invariant. 
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Here we have defined the magnetic helicity spectrum g(k,t ); the functions g±(k,t ) describe 
the relative weight of a single CK state W^ n ( x ; k ) with a definite helicity: 


s±(M) = Jh<4.(M)l 2 - 


(3.7) 


The energy of magnetic field can be related to g±(k ): 



Comparing (3.6) with (3.8), we find that the energy cost for a CK state W ± {x, k ) to carry 


one unit of helcity is k. 


B. Qualitative picture of the anomaly-driven inverse cascade of magnetic helicity 


We now ready to discuss the qualitative picture of the inverse cascade driven by anomaly. 
Let us first define the fermionic helicity: 



(3.9) 


where n A = j\ is the density of axial charge. From the anomaly equation (2.2) and the 


Maxwell equations, the total helicity of the system ho is conserved: 


(3.10) 


h 0 = h m + h F = const. 


Therefore, the system will tend to minimize the energy cost at fixed helicity. From the 


definition of “fermionic helicity” (3.9) and linearized equation of state n A = XI 1 A where x 


is the susceptibility, we observe that the energy per fermonic helicity is a a- On the other 
hand, as we discussed in Sec. |III A[ the energy per magnetic helicity for a single CK mode is 
k. Therefore for a hard (positive, i.e. of right circular polarization) CK mode k > a A (where 
without a loss of generality, we take a A to be positive as well), carrying helicity by chiral 
fermions is energetically favorable and consequently the helicity will be transferred from hard 
CK modes to chiral fermions. In contrast, the soft CK modes with k < a A are energetically 
favorable compared to chiral fermions, and fermionic helicity will be transferred to the soft 
components of magnetic helicity. Because the total helicity is conserved, this transfer will 








deplete the value of a A , and so the transfer of fermionic hclicity will occur to softer and 
softer CK modes - therefore, we find an inverse cascade of magnetic hclicity. 


We are now ready to extend our discussion of individual modes to the evolution of the 
entire hclicity spectrum g(k,t). For definiteness, let us assume that the total hclicity of the 
system h 0 is positive. It is convenient to introduce a characteristic energy scale kh associated 
with total hclicity: 


kh = 


C\h 0 


XV ’ 


(3.11) 


with V the volume. The quantity kh can be interpreted as the energy per hclicity if the entire 
hclicity h o is carried by chiral fermions. The fate of the system depends on the values of kh 
and k min , the lowest possible eigenvalue of a CK state allowed by the boundary conditions. 
If kh < fc min , eventually all magnetic hclicity will be transferred to fermonic helicity. In 
contrast, if kh > k m i n , the hclicity will eventually be carried by magnetic fields and the 
configuration of magnetic held will approach a single CK state H 7 ^ (fc m i n ) that minimizes 
the energy at a fixed helicity. We will confirm this scenario by a quantitative analysis in the 
next section. 


IV. MAXWELL-CHERN-SIMONS EQUATIONS 


The CME current can be described by adding the Chern-Simons term to the Maxwell 
theory [25] , Assuming that the gradients of chirality distribution are small, the resulting 
sourceless Maxwell-Chern-Simons (MCS) theory acquires the usual Maxwell form with an 
extra term in the current describing the CME: 


„ „ dE 

VxB = w 


V x E = 
V B = 0, 


+ Jem 

dB 


dt ' 

V -E = 0. 


Here, the electric current j*em includes the Ohmic and CME components: 

Jem = crE + a A B , 


(4.1a) 

(4.1b) 

(4.1c) 

(4.2) 


where a is the usual electrical conductivity. Let us now take the curl of the equation (4.1a) 


and use (4.1b) and (4.1c) to obtain 

ad t B(t , x) = X7 2 B + a a (V x B ) 


(4.3a) 
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In (4.3a), we have neglected d?B term - this should be a good approximation for time 
scales larger than 1/cr. Here, as mentioned above, we also neglect the spatial dependence of 
n A and relate n A to tia via the linearized equation of state ha = n a/x'i this yields 


&A(t) = 


CAn A (t) C'A f ,3 / ,N 

-— « — / d 6 xn A (x,t 

X XV J 


(4.3b) 


In accord with our assumption of small gradients of the axial density we will neglect the 
spatial component of axial current j A . The evolution of n A (t) is thus related to the evolution 
of Chern-Pontryagin density E ■ B by the anomaly equation (0: 


d t n A (t) = ~^f d 3xE ' B ■ 


(4.3c) 


Eqs. (4.3) give the simplified version of the MCS equations that we about to solve. 


A. General solutions 


Once we apply the decomposition (3.3), (4.3a) becomes a differential equation describing 
the time dependence of ai m (k,t ): 


S<“i™(M) = <r 1 [-k 2 ±a A (t)k}af m (k,t). 


(4.4) 


From the definition of chiral magnetic conductivity a A = C a ^a and (3.10), we have 

hm(f) 


0A(t) = h 


hr i 


(4.5) 


where kh is defined in (|3. 11). The solution to (4.4) and (4.5) can be obtained as follows 


(see also Refs. mm- First, integrating (]4.4[), we get 


a 


(M) = expl <r 1 [-kh±k0(t)} h 


(4.6) 


where af m r (k) = = 0, k) is determined by the initial value of magnetic field and 

9(t)= [ dt'a A (t'). (4.7) 


Now, from the definition (3.7) we have 


g ± {k, t) — gf{k) exp < 2cr 1 [~k 2 t ± k 9(t)] l, 


(4.8) 
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where gf(k) = gf(t = 0 ,k) denotes the initial magnetic helicity spectrum. Finally, 9{t) 


— „± 


(and thus a a (t)) will be determined from the consistency condition (4.5): 


1 r°° rlk 

a A (t) = k h { 1 - — / —k[g + (k,t) - g_(k,t)} 

"-o Jo ^ 


(4.9) 


Before presenting numerical solutions, we now discuss the evolution of individual CK 
modes ci^JJJk, t) as described by (4.4). Without losing generality, let us assume that a A (t) > 


0. Then the negative helicity mode a^ m (k,t) will decay exponentially, so let us concentrate 
on the evolution of the positive helicity mode a’i m (k,t). For hard modes, i.e. the modes 
with momenta k a A , a'i m (k,t ) decays exponentially exp(—cr _1 £; 2 t), as usual, ffowever, 
the soft helicity mode k < a A (t) will grow exponentially. This unstable mode has been 
noticed before in various contexts, see Refs. pH mmm for examples. The growth of 
soft CK modes could be anticipated from the discussion in Sec. |111B| the system tends to 
minimize the total energy while preserving the total helicity, and the soft CK state possesses 
the lowest energy at a fixed helicity. 


V. THE INVERSE CASCADE OF MAGNETIC HELICITY 


A. The initial conditions and Hopfion solutions 

As discussed in the previous section, the evolution of 9(t), a A (t) and g(k,t ) can be de¬ 
termined once the initial condition for the configuration of electromagnetic field is specified. 
Since we would like to investigate the evolution of topology of magnetic flux, we assume 
that initially the electromagnetic field with a non-zero magnetic helicity is localized at a 


short spatial scale much shorter than k, 1 defined in (|3.11|). We therefore take the Hopfion 


solution j22] to vacuum Maxwell equations as the initial configuration. This solution carries 


non-zero helicity (which we assume to be positive) that is due to the second term in (2.1) and 
a finite energy. It may be interpreted as a soliton wave solution to the Maxwell equations. 
For a Hopfion solution with a total initial magnetic helicity h m j, the electromagnetic held 
can be expressed in terms of CK states as [27] (see also Ref. 


B 


Hopf 


[X,t) = 


4 h 


m,I 


37T 


^HopfOM) 


4 h 


m.I 


37r 


dkk 2 e~ kLEM [(kLl M )W+(x,t)e~ ikt + c.c.] , 
dkk 2 e~ kLEM [(— ikLft M )W^_(x, t)e~ lkt + c.c.] , 


(5.1a) 

(5.1b) 
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where Lem characterizes the size of the Hopfion. Consequently, only modes af^kjt) are 
non-vanishing (see (3.3)) and the initial magnetic helicity spectrum is given by: 


9i(k) = 2kL ™ 


(5.2) 


The peak of gi(k) at k p = 1/Lem defines a characteristic size Lem of the configuration. 


B. Stages of the inverse cascade evolution 


We are now ready to solve the evolution equations (4.4) and (4.9). In addition to 


the initial condition discussed in Sec. V A the evolution also depends on the dimensionless 
ratios LEukh,h m j/ho,<r/kh . We would like to model the situation in which initially the 
helicity is dominated by the contribution from the EM field, thus h m ~ ho and Lem > k^ 1 . 
To be concrete, in this subsection we present the results corresponding to the solution 


with (L EM kh,h m j/h 0l <7/kh) = (1/2, 0.8, 0.4). We have also numerically solved (4.9) with 
different choices of Lem^, h mi //h 0 , cr/fc^; the results are qualitatively similar. 

In Fig. [2], we plot the time dependence of the chiral magnetic conductivity a Ait)-, the peak 
of g(k,t), Lpeak(t), and magnetic helicity h m (t). From Fig. [2j we observe that according to 
the behavior of cr^(t), /c pea k(t), h m (t), the evolution of system can be generally divided into 
the following three stages listed below. For reference, we also plot the magnetic helicity 
spectrum g(k,t ) at initial time t = 0, and three representative times corresponding to the 


three stages in Fig. 3(a) 


1. Stage I: in this stage, the magnetic helicity h m {t ) decays exponentially. Due to the 
conservation of total helicity, the magnetic helicity is transferred to fermionic helicity 
- thus we observe a fast growth of a Ait). Meanwhile, k peak {t) starts decreasing but 
is still larger than oy^f). Stage I ends when magnetic helicity h m (t) becomes small 
and ct.a(£) is close to k h . The duration of “stage I”, which we denote as 77 , can be 
estimated from the decay rate of magnetic helicity in this stage, cr _1 L| M , as indicated 


by (5.5). We therefore have: 


77 ~ <7 L 


EM 


(5.3) 


2. Stage II: in this stage, the total helicity ho is dominated by fermionic helicity hp. In 
other words, a Ait) approximately equals to k h and we observe from Fig. [2] that a Ait) 
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5.00 


Stage I Stage II Stage III 



FIG. 2. (Color online) The time dependence of chiral magnetic conductivity a A (t) (red), the peak 
of magnetic helcity spectrum k peak (t) (blue) and magnetic helicity h m (t ) (green). The schematic 
divisions of three stages (see text) for the evolution of the system are sketched in dotted horizontal 
lines. Dashed black curve illustrate f^ 1 / 2 asymptotic behavior of cr A (t), k peak (t). Black dotted curve 
below the green curve plots h m (t ) by solving Maxwell’s equation in the absence of anomaly effect. 


The numerical results are determined by solving (4.4) and (4.9) with (LvMkh,h m ,i/ho,cr/kh) = 
(1/2, 0.8, 0.4). 


changes slowly, while k peak (t) continues to decrease. “Stage II' 7 ends when K pea k(£) is 
close to a a- 


3. Stage III: in this stage, both a A (t) and k peak (t) decrease. The fermionic helicity is 
transferred to magnetic helicity and eventually h m (t ) will approach ho. At late times, 
cr A (t) ~ k peak (t). This corresponds to the configuration in which the energy cost per 
helicity for fermionic helicity is approximately equal to that of magnetic helicity. In 
this case, the following relation holds V x B « k peak (t)B ps a A (t)B. It is clear from 
the log-log plot Fig. [2] that cr^(i), k peak (t) behave as a power law in time t: 

k h (t ) « a A it) oc t~ p . (5.4) 

Meanwhile, the evolution of g(k,t ) becomes self-similar: 

g(k,t) ~t a g(t fi k), (5.5) 
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(a) (b) 

FIG. 3. (Color online) The evolution of magnetic helicity spectrum g(k,t). (Left): g(k,t) at 
initial time t = 0 (red) and three representative time (corresponding to stage I, II, III) respectively. 
(Right): t~ a g(k,t ) vs t^k in the self-similar stage (stage III) of the evolution. 


where g (t^k) is the scaling function and 


a = 1, 0 = 1/2, 


(5.6) 


are scaling exponents. It is easy to see that once g(k,t) becomes self-similar as in 


(5.5), fcpeak(^) is determined by the peak of the scaling function g(t^k). Therefore (3 
in (|5.5l) is identical to (3 in (5.4). 


In Fig. 3(b) , we plot g(k,t)/t a vs tPk with scaling exponents given by ( |5.6[ ) for 


different t in Stage III. The self-similar behavior of g(k,t) is evident from Fig. 3(b) 


At this point, critical exponents (5.6) are found numerically. In Sec. VC, we will 


determine the scaling function g and derive (5.6) analytically. 


To close this section, we emphasize that chiral anomaly plays a crucial role during the evo¬ 
lution of chiral magnetic conductivity cr^{t) and of the magnetic helicity spectrum g(k,t). 


Indeed, without the CME current term in Maxwell equation (4.1a) and with no trans¬ 
fer of helicity between magnetic held and chiral fermions, g(k,t ) would simple decay as 
exp(— 2cr~ 1 k 2 t) (see also black dotted curve in Fig. [ 2 ]) and self-similar behavior would be 
absent. 
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C. Self-similar evolution and scaling behavior of g(k,t) 


In this section, we would like to understand the origin of the scaling exponents (5.6) 
found numerically, and to determine the scaling function g(t^k) . 

First of all, we note that a, (3 are not independent. As h m & h 0 at late time (c.f. Fig. [2j) , 
we have: 

h m {t) = I ^-kg(k) = J ^-kt a g(t p k) = t a ~ 2/3 J ^-g{x) ~ h 0 , (5.7) 

where we have introduced a new variable: x = t^k. We therefore have: 


a = 2/3. 


Moreover, if (4.8) can be matched to the scaling form (5.5), we must have 


> 9 == 


Consequently, for self-similar evolution, we have from (4.7): 


9(t) = 2<TA(t)t . 


Substituting (5.10) into (4.8), we obtain: 


g(k , t) = gi(k ) exp 2 a 1 [- k 2 + 2 a A (t)] t 


= gi(k ) exp 


2a\(t)t 


a 


exp < — 2a 1 [k — cr^f)] t > . 


(5.8) 


(5.9) 


(5.10) 


(5.11) 


If the width of the Gaussian in (5.11 ) is sufficiently narrow, we further have /c pe ak(^) ~ cta(^) 
and 

g(k,t) oc 5 , /(fcpeak(^))exp | - 2cr -1 [k - A: P eak(t)] 2 t| • (5.12) 

To summarize, the system will spend a long time at the stage of self-similar evolution. 
In this stage, k pea k{t) decreases as t -1 / 2 . This implies that a large-scale helical magnetic 
held will develop. With the growth of t, the width of the Gaussian becomes more and more 
narrow, and g{k,t) will become proportional to delta-function : 


g(k, t) ->■ S[k- k peak (t)} . (5.13) 

In this limit, the system is described by a single CK state IFp (fc pe a,k(t), t). Eventually, the 
evolution will end when /c pea k(^) ~ l/L where L is the size of the system. Here we have found 
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self-similar evolution by solving the MCS equation with the Hophon initial condition (5.1). 
However, as our analysis does not rely on any particular feature of the Hophon solution, we 
expect that self-similar evolution is a general feature that at late times does not depend on 
the choice of initial conditions. 


VI. IMPLICATIONS FOR HEAVY-ION COLLISIONS 


Let us now establish whether a self-similar evolution of magnetic helicity can be realized 
in experiment; here we will focus on QCD matter created in heavy-ion collisions. As we 


discussed in Sec. Ill B to realize the inverse cascade of magnetic helicity, the energy cost 


per helicity if total helicity is carried by chiral fermions kh (see (3.11)), should be larger 
than fc min , the minimum eigenvalue of a CK state. As k m ; n is of the order 1/L, where L 
is the size of the system, we need to check whether kh > 1/L with L ~ 10 fm can be 
realized in a heavy-ion collision. As both magnetic helicity and fermionic helicity would 
contribute to the total helicity, we will estimate their contributions separately. In heavy-ion 
collisions, the initial axial charge density can be generated by sphaleron transitions and/or 
by the color flux tubes during the early moments of the heavy-ion collision |29j . To make a 
rough estimate, we will follow Ref. ra and consider an (optimistically large) value of the 
chiral chemical potential /i^, of the order of 1 GeV. The resulting a a is then of the order 
0.01 GeV. If the total helicity ho originates mostly from this initial axial charge, 1 /kh would 
be at least of order 20 fm, which is much larger than the typical size of the fireball L ~ 10 fm 
created in a heavy-ion collision. This estimate is also in agreement with Ref. [3D] in which 
the relevance of CK state to heavy-ion collisions was discussed. To summarize, in order to 
satisfy kh > 1/L, initially the dominant contribution to the total helicity should be from 
magnetic helcity h m . 


We now estimate magnetic helicity in a heavy-ion experiment. We first note that from 
the EM field pattern created by spectators in heavy-ion collisions (c.f. Fig g ), e b is 
positive in the upper half region, and negative in the lower half region. The EM field would 
thus be helical, with opposite helieities in the upper and lower half region. To estimate the 
magnitude of this magnetic helicity, we assume \B\ & \E\ and in RHIC A ~ E Turkic 
where tb,rhic is the typical lifetime of magnetic field at RHIC. The typical (peak) strength 
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FIG. 4. (Color online) A typical configuration of EM field in the transverse plane in non-central 
collisions. The two circles indicate the edge of the colliding nuclei. The solid line (red) shows the 
magnetic field B and the dashed lines (blue) show the electric field E. The inner product between 
electric and magnetic fields E ■ B becomes positive (negative) on the upper (lower) side of reaction 
plane (RP). 


of magnetic field at RHIC is 

cRrhic = cp/rn 2 ^ , (6.1) 

where Cp varies in the range of 1 to 10 depending on the impact parameter after event 
averaging, but in a given event can be significantly larger than this average value due to 
fluctuations. We therefore have: 

C'a f d 3 xA RHW • -Brhic e ~ 2 C\ (ei?RHic) 2 Tp ,RHIC 


k 


ft,RHIC 


XaV 


N fXf 


= OEM (Cem) ( 


N 2 \ ( eB RHIC ) Tp tRH ic 


= 1.0 x 10 b c 2 pTp tRHI c fm 2 . (6.2) 


N f Xf 


Here we have assumed that the axial susceptibility is xa = NfXf where X/ is the quark 
number susceptibility known from the lattice measurements PH |22] to be Xf ~ 1-0 T 2 for 
temperatures higher than T c . We consider the case of u, d flavors contributing to the CME 
current and thus put Nf = 2, Cem = 5/9. We also assume that at the initial stage of 
heavy-ion collisions T 0.4 GeV. To compute kh at LHC, we further take: 


ei?LHC 

ei^RHic 


7lhc 
7 RHIC 


13 g Tp, LHC _ f 7LHC \ 
TB,RHIC \7 RHIC J 


-1 


(6.3) 
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We then get: 

&/t,LHC ~ ( ^ LHC ^ /ch.RHic = 1-4 x 10~ a c 2 b tb,kric & n 2 • (6.4) 

\7rhic / 

For the purpose of estimate, we will take tb, rhic = 1 f m - Therefore in order to satisfy 
kh > 1 /L ~ 0.1 fm -1 , we need to select events with cb ~ 100 at RHIC and cb ~ 26 at LHC. 

We conclude that observing the self-similar cascade of magnetic helicity in heavy ion 
collisions will be challenging. However the estimated magnetic helicity is quite large, and is 
likely to affect the evolution of the quark-gluon plasma. It can lead to interesting observ¬ 
able effects. For example, since the sign of magnetic helicity is different in the upper and 
lower hemispheres (above and below the reaction plane), the decay of the magnetic held at 
freeze-out will yield the photons with opposite circular polarizations. Since the direction of 
magnetic held, and thus the signs of magnetic helicity, can be determined experimentally 
by measuring the spectators, one can measure the polarizations of photons by summing 
over many events (the sign of polarization will not huctuate event-by-event). The pho¬ 
ton polarization can be measured through photon conversion into e + e _ pairs by extracting 
the angular distribution of the electrons and positrons. We believe that the observation of 
these opposite circular polarizations of the produced photons will be a unique signature of 
magnetic helicity in heavy ion collisions. 


VII. SUMMARY 


To summarize, the chiral anomaly couples the evolution of axial charge density and 
electric-magnetic (EM) held in the plasmas possessing chiral fermions. By solving the 
Maxwell-Chern-Simons equation in the presence of CME current, we analyzed the real time 
evolution of the magnetic helicity spectrum. We initialized the system by assuming that it 
contains a seed of a helical magnetic held at short distances, with helicity carried by the 


linked magnetic held configuration, i.e. by the second term in (2.1). As summarized in 


Fig. [2j we found that the magnetic helicity hrst gets transferred to fermionic helicity and 
then fermionic helicity is transformed back into magnetic helicity, but this time to self-linked 


Chandrasekhar-Kendall (CK) configurations characterized by the second term in (2.1). We 


have argued that the CK states that minimize magnetic energy at a hxed helicity represent 
the final stage of the magnetic helicity evolution. We found that at late stage, this evolution 
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becomes self-similar, and describes the growth of the self-linked CK state. 


The role of fermions is to mediate the magnetic reconnections that are necessary to 
transfer the magnetic helicity from the linked to self-linked configurations of magnetic flux, 


see (2.1). We expect that our hirelings apply to all systems that possess the CME current. 
In addition to the quark-gluon plasma discussed above, the growth of magnetic helicity can 
be expected in Dirac semimetals that exhibit the CME in parallel electric and magnetic 
fields [33]. Experimentally, this generation of magnetic helicity can manifest itself through 
the emission of circularly polarized photons in the THz frequency range characteristic for 
Dirac semimetals 


As a natural extension of this work, it will be interesting to study the inverse cascade in 
the framework of anomalous MagnetoHydroDynamics (MHD). In this case, the anomaly can 
couple the kinetic helicity carried by the fluid, magnetic helicity and fermionic helicity. While 
the inverse cascade of magnetic helicity is a traditional topic of magnetohydrodynamics, the 
role played by the chiral anomaly has not yet been fully explored. As another extension of 
this work, it would be interesting to include the spatial dependence of axial charge density 
in the MCS equations. This would allow us to study the evolution of domains with non-zero 
axial charge throughout the inverse cascade of magnetic helicity. 
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Appendix A: A useful representation of Chardrasekhar-Kendall states 


A single CK state W ± (x] k) can be represented as the linear combination of toroidal field 
T(x;k) and poloidal field P(x;k) 

W± n (x; k ) = Ti m {x ; k) =F iPi m (x; k ), (Al) 

where 

Tim{x ; k) = ji{kr)Xi m {9 , 0), Pi m (x; k) = f V x T] m (cc; /c), (A2) 

k 

where ji{kr) denotes spherical Bessel functions. Here, we have defined: 

Xim(0, 0) = 0)] , i = -i (rxV). (A3) 

In the above equations, Yi m (9, 0) is the usual speherical harmonic functions. 
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